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Abstract — Consider a multiuser system where an arbitrary 
number of users communicate with a distributed receive array 
over independent Rayleigh fading paths. The receive array 
performs minimum mean squared error (MMSE) or zero forcing 
(ZF) combining and perfect channel state information is assumed 
at the receiver. This scenario is well-known and exact analysis is 
possible when the receive antennas are located in a single array. 
However, when the antennas are distributed, the individual links 
all have different average signal to noise ratio (SNRs) and this 
is a much more challenging problem. In this paper, we provide 
approximate distributions for the output SNR of a ZF receiver 
and the output signal to interference plus noise ratio (SINR) of 
an MMSE receiver. In addition, simple high SNR approximations 
are provided for the symbol error rate (SER) of both receivers 
assuming Af-PSK or M-QAM modulations. These high SNR 
results provide array gain and diversity gain information as well 
as a remarkably simple functional link between performance 
and the link powers. These metrics can be used by a scheduler 
tasked with the job of selecting a subset of users for simultaneous 
reception, when only long term channel state information is 
available. 

Index Terms — Macrodiversity, MMSE, ZF, Outage probability, 
Optimum combining, Zero-Forcing, Network MIMO, CoMP. 

I. Introduction 

With the advent of space diversity systems, decoupling 
users through channel aware signal processing techniques in 
the presence of multiple access interference (MAI) and noise 
has become an integral part of the system design. There are 
various processing techniques now widely adopted in research 
and standards Among them, linear combining methods 
are popular for their simplicity despite the fact that they are 
not optimum in a maximum likelihood sense. Two key linear 
combiners are zero forcing (ZF) and minimum-mean squared- 
error (MMSE). Although they are not optimal, the MMSE 
receiver satisfies an alternative criterion, i.e., it minimizes the 
mean squared error (MSE) and ZF is known to eliminate MAI 
completely. 

The performance analysis of such linear receivers is of great 
interest in wireless communication (2) as it provides a baseline 
link level performance metric for the system. Today, perfor- 
mance results for MMSE/ZF receivers are well known for 
microdiversity systems where co-located diversity antennas 

D. A. Basnayaka, P. J. Smith and P. A. Martin are with the Department of 
Electrical and Computer Engineering, University of Canterbury, Christchurch, 
New Zealand. E-mail:{dush, p. smith, p.martin}@elec.canterbury.ac.nz. 

D. A. Basnayaka is supported by a University of Canterbury International 
Doctoral Scholarship. 



at the base station communicate with distributed users 0, 
lIU, 0. Macro-scale diversity combining has recently become 
more common from a variety of perspectives |6|, 0. Any 
system where both transmit and receive antennas are widely 
separated can be interpreted as a macrodiversity multiple 
input multiple output (MIMO) system. They occur naturally 
in network MIMO systems 0, and collaborative MIMO 
concepts p. 69] and iflOl . 

The performance of macrodiversity systems has been inves- 
tigated via simulation fTD . but very few analytical results 
appear to be available. The reason for the lack of results is the 
complexity of the channel matrix that arises in macrodiversity 
systems. Some analytical results are available for the dual user 
case in lTT2l for macrodiversity MMSE and ZF receivers. In 
[ 12 1, they consider the statistical properties of the output signal 
to interference plus noise ratio (SINR)/signal to noise ratio 
(SNR) of MMSE and ZF receivers respectively and obtain 
high SNR approximations of the symbol error rate (SER). 
In IPT31 . the SER performance of macrodiversity maximal 
ratio combining (MRC) has been exactly derived for arbitrary 
numbers of users and antenna configurations. Ergodic sum 
capacity of macrodiversity MIMO multiple access channel is 
considered in [14] and derive very tight approximation for 
ergodic sum capacity in a highly compact form. Rayleigh 
fading is assumed in fl2l . |[T3l and lfl4ll . 
In this paper, we extend the results in Ifl2ll to more general 
user and antenna configurations. In particular, the contributions 
made are as follows: 

1. We derive the approximate probability distribution func- 
tion (pdf) and cumulative distribution function (cdf) of 
the output SINR/SNR of MMSE/ZF receivers. 

2. The approximate cumulative distribution functions are 
shown to have a remarkably simple form as a generalized 
mixture of exponentials. 

3. High SNR approximations for the SER of MMSE/ZF 
receivers are derived for a range of modulations. 

4. The SER results are used to derive diversity order and 
array gain results. 

5. The high SNR results are simple, have a compact form 
and can be used to gain further insights into the effects 
of channel distribution information (CDI) on the perfor- 
mance of macrodiversity MIMO systems. 

The rest of the paper is laid out as follows. Sec. II describes the 
system model and receiver types. Sec. Ill provides preliminary 
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results which will be used throughout the paper. The main 
analysis is given in Sees. IV and V. Sees. VI and VII give 
numerical results and conclusions, respectively. 

II. System Model 

The multiuser MIMO system investigated in this paper 
consists of N distributed single antenna users communicating 
with rift distributed receive antennas in an independent flat 
Rayleigh fading environment. The C nR 
given by 



xi 



receive vector is 



Hs 



n. 



(1) 



where the C Nxl data vector, s — (si, s 2 , ■ ■ ■ , sn) T , con- 
tains the transmitted symbols from the N users and it is 
normalized, so that E'-jJsil 2 } = 1 for i = 1,2,. ..,N. 
n is the C nRXl additive-white-Gaussian-noise (AWGN) vec- 
tor, n ~ CM (0, cr 2 /), which has independent entries with 
^{l n j| 2 } = °' 2 ' f° r * = 1> 2) . . • , TiR. The channel matrix 
contains independent elements, Hi k ~ CA/"(0, Pik), where 
E{\H ik \ 2 } = P ik . A typical macrodiversity MU-MIMO 
multiple access channel (MAC) is shown in Fig. Q] where it is 
clear that the geographical spread of users and antennas creates 
a channel matrix H, which has independent entries with 
different P ik values. We define the C UrxN matrix, P = {P ik }, 
which holds the average link powers due to shadowing, path 
fading, etc. 

By assuming that perfect channel state information is available 




Fig. 1. System diagram. To reduce clutter, only paths from a single source 
are shown. 

at the receiver side, we consider a system where channel 
adaptive linear combining is performed at the receiver to 
suppress multiple access interference (TJ. Therefore, the C Nxl 
combiner output vector is f = V H r, where V is an C nRXN 
weight matrix. In this work, we consider two well-known 
linear combining schemes: MMSE and ZF. The structure 
of V and the resulting output SINR/SNR for MMSE/ZF 
schemes are well-known and are given below. Without loss 
of generality, we assume that the index of the desired user is 
i = 1, The combining vector and output SINR of the MMSE 
receiver for user 1 are given by 0, lfT31 as 



R hi , 



(2) 



SlNR = h"R~ 1 h 1 , 



(3) 



where 



JV 



R = Y,hkh^ + a 2 I, (4) 

where H = (h\,h2, ■■■ , ft/v), Defining «2, ... ,Vn similarly 
gives V = (t>i,«2, • • • ,«a'). The vectors, h k , clearly play 
an important role in MMSE combining and it is useful to 
define the covariance matrix of h k by P k = E ^h k h k ^ = 

diag(Pi fc , P 2k , P nR k)- From 0], 0, the combining ma- 
trix, V, and output SNR of the ZF receiver for tir > N are 
given by 

(5) 



H ( H H H 



and 



SNR = 



1 


a 2 


\h b h)~ v 








li 



(6) 



where indicates the (1, l) th element of matrix B. 



III. Preliminaries 

In this section, we derive some useful results which will be 
used extensively throughout the paper. 

Lemma 1. M4V Let X be an n x n random matrix with, 
(E{\X^) E{\X 12 \ 2 } ...\ 
E{\X 21 \ 2 } E{\X 22 \ 2 } ... 



A = E{XoX} 4 



Vi^l^l 2 } E{\X n2 \ 2 } ...J 



(7) 



where o represents the Hadamard product. With this notation, 
the following identity holds. 



E { X"X } = perm (A) , 



(8) 



where perm{.) is the permanent of a square matrix defined in 

Proof: From the definition of the determinant of a generic 
matrix, X = {X i<k } i k=1 n , we have 



E- 



{\x H x\} 



(7 i—1 

n 



i=i 



(9) 



where a — (a\, a 2 , . . . , a n ) is a permutation of the integers 
1, . . . , n, the sum is over all permutations, and sgn (a) denotes 
the sign of the permutation. Since all the elements of X are 
independent, the only terms giving non-zero value expectations 
are Yi7=i X ai ,i Ti7=i Xbt.i, where permutation a = b. Hence, 



3 



using the permanent definition in |[T6l we have 

n 

=J2Ii A ^ -Penn(i4). (10) 



Corollary 1. [14] Let X be an m x n random matrix with, 
E {X o X} = A, where A is an m x n deterministic matrix 
and to > n. Then, the following identity holds. 



E- 



{ X"X } = Perm {A) 



(11) 



where Perm (A) is the permanent of the rectangular matrix A 
as defined in H16V . 



Proof: Using the Cauchy-Binet formula for the determi- 
nant of the product of two rectangular matrices, we can expand 
X H X as a sum of products of two square matrices. Each 
product of square matrices can be evaluated using Lemma Q] 
The resultant expression is seen to be the permanent of the 
rectangular matrix A, which completes the proof. ■ 

Corollary 2. M4V Let X be an m x n random matrix with, 
E {X o X} = A, where A is an m x n deterministic matrix 
and 777 > n. If the m x m deterministic matrix S is diagonal, 
then the following identity holds, 



E- 



| X H Y,X | = Perm (SA) 



(12) 



Proof: The result follows directly from Lemma [T] and the 
fact that E 5 o S 3 = £ for any diagonal matrix. ■ 
Next we give a definition for the elementary symmetric 
function (esf) of degree A: in 77 variables X\ , X2 , ■ ■ ■ , X n from 
fiin . Let e k [X U X 2 , . . . ,X n ) be the k th degree esf, 



e-k \Xi,X2i ...,X Tl 



E 

l<h<l 2 <-<lk<n 



X, 



(13) 



It is apparent from (foi l that eo (Xi, X2, • • • , X n ) = 1 and 
e„ (Xi,X 2 , . . . , X n ) = XiX 2 ■ . ■ X n . In general, the esf of 
degree k in n variables for any k < 77 is formed by adding 
together all distinct products of k distinct variables. 

Lemma 2. H18V let X be an n x n complex symmetric 
positive definite matrix with eigenvalues Ai, . . . , A„. Then, the 
following identity holds. 



efc(Ai,A 2 ,...,A„) =Tr fe (X) , 



where 



Tr fc (X) 




1 < k < 77 
k = 
fc > 77, 



(14) 



(15) 



where <7 k. n is an ordered subset of {77} = {1, . . . , 77} of length 
k and the summation over all such subsets. X ae n denotes the 
principal submatrix of X formed by taking only the rows and 
columns indexed by ai :Tl . 

In general, X%.*' n denotes the submatrix of X formed by 
taking only the rows and columns indexed by <t^„ and „ 



respectively, where a^ n and \i^ n are length £ subsets of 
{1,2,..., 77}. If either ag n or [ig n contain the complete set, 
the corresponding subscript/superscript may be dropped. When 
&t,n — Hi,™ on ly one subscript/superscript may be shown for 
brevity. 

Next, we present three axiomatic identities for permanents 

ma. 

• Axiom 1: Let A be an arbitrary 777 x n matrix, then 

£ Perm ((A)' -) = $>erm ((A) ctq J = 1. (16) 

• Axiom 2: Let A be an arbitrary m x n matrix, then 

£Perm((A) CTtm ) = £ Perm ((A)^") . (17) 

a a 

> Axiom 3: For an empty matrix, A, 

PermL4) = l. (18) 



IV. ZF Analysis 



In this section, we derive an approximate cdf for the output 
SNR of a ZF receiver, a high SNR approximation to SER 
and also consider some special cases. The following pdfs for 
the columns of the channel matrix are used throughout the 
analysis. 



f(hk) 
for k = 1,2,...,N. 



1 

TT nR LP* I 



(19) 



A. CDF Approximations 

The output SNR of a ZF receiver in (|6]l can be written as 
Z = \hf (i - H 2 (H^H 2 ) _1 hA fci (20) 



where M 



G 

I H 2 (H^H 2 



(21) 



Hn and H 



Following the analysis in [12|, the characteristic function (cf) 
of Z is given by 



0z(*) 



E{e^}=E^ h " Mh >y (22) 

Conditioning on H 2 , the expectation over h\ in ( f22l can be 
solved as in lfl2l to obtain 

- 1/ Amp A' (23) 
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Substituting for M in ((23) gives 

^(<|H 2 ) - 



1 



1 



IDI 





H i 2 i H 2 




\D\ 


H?H 2 


+ ^H^P 1 D~ 1 H 2 



(24) 



(25) 



(26) 



where D = I — -\jtP±. Simplifying 

<P z {t\H 2 



gives 





H 2 H 2 




\D\ 


H%D~ X H 2 



(27) 



The full cf can then be obtained by averaging the conditional 
cf in d27i l. to give 



\D\ 



H 2 H2 



(28) 



where expectation is over H 2 . An exact analysis of (128l l is 
extremely cumbersome. However, for the dual source scenario 
where N — 2, d28l l can be solved in closed form [22|. Even for 
N = 2, the resulting exact expressions are complex. Hence, 
for arbitrary N we use a Laplace type approximation as in 
lfT2l to approximate and simplify the cf. This approximation 
has some motivation in the work of l23l and lfl4l . 11241 . It can 
also be thought of as a first order delta expansion [25 1. This 
approach gives 



<f>z(t) 



Applying Lemma Q] to ( |29l ) gives 

Perm(Q 2 ) 



1 H 


H 2 H 2 


} 


\d\e{ 


H*D l H 2 


} 



<Pz(t) 



(29) 



(30) 



|D|Perm(D -1 Q 2 

where the njj x (JV — 1) matrix Q 2 is defined by Q 2 = 
E {H 2 o H 2 }. P = (pjQ 2 ), and Pl = (P u , P 21 , . . . , P nRl f '. 
From Appendix |A] the denominator in ( f30b can be expanded 

as 

L 

\D\ Perm (l^Q,) = ]T (-jt)* <p f , (31) 



i=0 



where 

(pi = 



^l) 



perm 



((<9 



{N-l} 
Zlcrjsr-x.n, 



(32) 



Since perm ((Q^^ ^ ^ is independent of it is clear from 
that |D|Perm(D _1 Q 2 ) is a polynomial in t of degree 



L. Hence, (1301 becomes 



Perm(Q 2 ) 
W) w 7^1 



Perm(Q 2 ) 

Etc (&) ("J*)' 

Perm (Q 2 ) 

^ nf=i (w< - jt) 

Perm(Q 2 ) ^ ^ 



(33) 
(34) 

(35) 

(36) 



where tDj > for all i from Descartes rule of signs for the 
polynomial in fl34"l > and 



Vi = 



(37) 



It is clear from ( f30b that 0jj(O) = 1, since D = I when t = 0. 
Therefore, the cf will produce a valid pdf after inversion 1191 . 
From 1 20 1, the pdf and cdf of Z are given by 

1 



(38) 



-00 

i>Z p 00 



F z( z ) = 7T I / ^z(t>- Jtx dtdx. (39) 
Substituting (f36t in ( f38T > gives the approximate pdf 



Perm (Q 2 

2TT!fi L 



L 



-dt. (40) 



~ Jt 

Applying the following integral identity from lETTl . 

p > 
p<0, 

[i?e(«) > 0, R&{fi) > 0] , (41) 



-dx = I T ^ 
10 



we obtain the approximate pdf of Z as 
and the approximate cdf of Z becomes 



Perm(Q 2 ) ^ 77, 



(42) 



(43) 



The final pdf approximation in d42l has a remarkably simple 
form as a generalized mixture of L exponentials where L = 
n R -N + l. 



B. High SNR Approximations 



The cf in d28l l is a ratio of determinants, where D = I — 
-^jtP\. As the SNR grows, a 2 —> and keeping only the 
dominant power of a 1 in (l28T l gives 

2 \ n R -N+l 



4>z(t)=K 



-Jt 



(44) 
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where 



K 



\Pi\ 



tE ■ 



H2H2 



H 2 Pi 1 H2 



(45) 



Following the mgf based approach in B6I . the SER of a 
macrodiversity ZF receiver can be evaluated for A/-PSK 
modulation as 



PS = - I My I - 



9 



sm z 9 



d6. 



(46) 



where A^^(s) = ^^(—js), g = sin 2 (tt/M) and T = 
( - M ~ / 1 ' >7r . Note that linear combinations of equations of the 
form given in d46b also give SERs for il/-QAM in the usual 
way [26 1. Substituting ( 144b in d46l > gives 



P: 



(Gai) Gd + o(~r dd ), 



(47) 



where o(.) is the standard little-o notation and the average 
SNR is 7 = 4y. The diversity gain and array gain in (l4Tb are 
given by 

/ „ _\ -1/(«R-JV+1) 

G a = (KqI) 



G d = n R - N + 1 
where I is given by 



sin 2 6 
9 



(n R -N+l) 



dO. 



(48) 



The high SNR expression derived in (l4Tb will be exact, if 
and only if Ka is exact. An exact calculation of Ka for the 
N — 2 case is presented in 11221 and shown to have a complex 
expression. This work suggests that in the general case an 
exact calculation is likely to be either excessively complicated 
or intractable. Hence, in this work, we use a Laplace type 
approximation for Kq in d45l > to obtain a more compact and 
insightful expression. Hence, we use the following approxi- 
mation 



K 



Using Lemma Q] d49l is given by 

Perm (Q 2 ) 



1 


h 2 h 2 


} 




H2 Pi 1 H2 


} 



K a 



\Pi\ Perm IP^Q 



(49) 



(50) 



Note that when N = 2, approximate Ka has simpler expres- 
sion |[T2l . which gives 



Tr(P 2 



\Pi\Tr[Pi 1 P 2 



(51) 



V. MMSE Analysis 

A. CDF Approximations 

In this section, we derive the approximate cdf of the output 
SINR of an MMSE receiver and a high SNR approximation to 
the SER. Let Z be the output SINR of an MMSE receiver given 
by (01. Following the same procedure as in the ZF analysis, 



the cf of Z is 

fe(t) = E{e 3tz ) =E^e jth " R ~ lhl Y (52) 



Next, the cf conditioned on H 2 becomes ifTZl 

1 



K*l#2) = 



I-jtR^Pi 



(53) 



Since R = a 2 1 + H2H2 , the conditional cf in ( 1531 becomes 



<pz(t\H 2 ) 



I-jt(<r*I + H 2 H%) Pi 





a 2 I + H 2 H% 






-jtPi+H 2 H% 



(54) 



(55) 



1 + XX 



H 



I + X H X 



Using the determinant identity, 
where the rank of the identity matrix is obvious from the 
context, in ( TSBl along with some simple algebra, we get 



*>z(t\H 2 ) = 





a 2 1 + H%H 2 




\D\ 


<j 2 I + H 2 H D 1 H 2 



(56) 



where D = I \]tPi- Note the similarity of (l56l with 

ITLZl eq. 14]. Then, the full cf can be solved by averaging the 
conditional cf in (1561 over #2. Hence, 



a 2 I + H^H 2 



a 2 I + H?D~ 1 H 2 



(57) 



Using a similar approach as in the ZF analysis we approximate 
d57l > to get 



1 E { 


a 2 I + H^H 2 


} 


P\e{ 


<j 2 i + iffrr 1 2J 2 


} 



(58) 



In Appendix [B] we obtain the expectation in the numerator of 



as 



N-l 



e{\* 2 i+ h^h 2 I }=EE Perm ((^r*"- 1 ) (^ 2 



fc=0 <T 



Appendix [C] gives the denominator of ( 1581 as 

\D\ E { \a 2 I + H H 2 D X H 2 } = £ (-jtf Vi 

8=0 

where 

AT-l 

¥>i = E ^ ife ( CT2 



. JV- 



-fc-1 



Af-fc-1 



(59) 



(60) 



(61) 
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Substituting (|59]l and (|60]l in we get, 

Ef=o 1 E.Penn((0 2 )^- 1 )(a 2 ) J 



where 



Et R o(-J*)Vi 
e(Q 2 ) 



^ fl Er=ote)(-i*r 

e (Q 2 ) 

nr=i - ' 

where > for all i from Descartes rule of signs and 

JV-l 

0(Q 2 ) = EE Perm W2)^ N - i )(^) 

fc=0 a 

The final expression for <fiz(t) then becomes 

Vi 



N-k-l 



where 



Vi T~I n R I 

As in the ZF analysis, the pdf and cdf of Z are given by 
1 



(62) 
(63) 

(64) 
(65) 

(66) 
(67) 



M*0 = ^ / Mt)e~ jtz dt 



Fz{z) 



— OO 

z poo 



(68) 

<pz{t)e' jtx dtdx. (69) 



2tt 7o i-c 

Using the identity in (PiTT i to solve d68l l and (|69l we get the 
approximate pdf of Z as 

©(Q 2 ) 



(70) 



and the cdf of Z becomes 



F z (z) = ^9Af^Hi(l-e-^ ). (71) 



In contrast to (142 b . where the ZF SNR is a generalized mixture 
of L exponentials, (f70l can be identified as a generalized 
mixture of > L exponentials. Since the MMSE SINR 
has more mixing parameters {nji rather than L) it might be 
expected that these increased degrees of freedom will result 
in a better approximation. Alternatively, the more concise ZF 
result may provide a simpler expression for use in system 
design and understanding. These issues are considered in Sec. 

ivm 



B. High SNR Approximations 



The cf in d58t is a ratio of determinants. As the SNR grows, 
a 2 — > and keeping only the dominant power of a 2 in 
gives 

2 \ n R -N+l 



4>(t) = K (-jt) 



-jt 



(72) 



Ko(s) 



\Pi\ 





H 2 ' H2 




H~2j 


Pi 1 H 2 


+ sl 



(73) 



Hence, from d46i >. the SER at high SNR becomes 

) K °i-^) de - (74) 

Wo \ 9 J \sm-ej 
As in the ZF analysis, an exact calculation of Kq appears 
difficult and we use the Laplace-type approximation again to 
give 



K (s) 



From Lemma Q] and j59) , we have 
Perm(Q 2 ) 



1 *{ 


H 2 H% 


} 




si + H^P^H 2 


} 



K (s) 



,N-i-l 



(75) 

(76) 

(77) 
(78) 



where Q — Ecr P erm ((^1 ^Qt^ ) and — §i are the 

roots of Eilo^i s"^ 1 ' Since CiV-i = Perm (p^Q^^j, 
Kq (s) in ( f78T > becomes 



Perm(Q 2 ) x - X» 

|Pi|Cjv-i ^^i + s' 



K (s) 



Perm(Q 2 ) 



7V-1 



where 



|P 1 |Perm(p^ 1 Q 2 ) 

(n-i 



From d74l > and f79) . we obtain 

P s °° = {G a *f)- G * + o(T Gi ) 



(79) 



(80) 



(81) 



where the diversity order and array gain, Gd and G a respec- 
tively, are given by 



and 



where 



G a 



G d = n a -N + 1, 

_Perm(Q 2 ) 
Pen 



IPilPermfP^ 1 ^ 

i(p) = — / ^ 



j(p) 



<" ^sin 2 0) i+1 



(82) 



(83) 



(84) 
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The integrals in (|84] > can be solved in closed form as in [22|. 
Hence, the final result becomes 



1 N ~ x / 



where 



c -2 m . 



1 f c sin' 



t Jo a + sin 



(85) 



(86) 



VI. Special Cases 

In this section we present the special case where = N, 
i.e., the system schedules as many simultaneous users as the 
number of receive antennas. In this particular scenario, the ZF 
cdf analysis in Sec. IIV-AI has an intriguing form. From d30l >. 
the cf of Z becomes 

Perm(Q 2 ) 



(87) 



\D\Perm[D^ 1 Q 2 
When riR — N, the denominator of d87t simplifies to give, 

\D\ Perm (£r x Q 2 ) = ^ (l - ^P a ) perm (Q i2 ) , (88) 

i— 1 ^ ' 

where Q i2 is Q 2 with the i* row removed. Then, (l87T i 
simplifies to 



Perm (Q 2 



Er=i Perm (Q l2 ) - f El=i PuP^m (Q, 
Perm(Q 2 ) 



(89) 



(90) 



Perm(Q 2 ) - ^perm(P)' 

Inverting the cf expression in ( |90l gives the approximate pdf 
of Z as the simple exponential 

j~ z {z)=oHe-° H \ (91) 

where 9 = Perm (Q 2 ) /perm (P). 

VII. Simulations and Numerical Results 

In this section, we simulate the macrodiversity system 
shown in Fig. [2] where three base stations (BSs) collaborate 
via a central backhaul processing (BPU) in the shaded three 
sector cluster. This simulation environment was also used in 
[12 1 and is sometimes referred to as an edge-excited cell. We 
consider the three BS scenario having either a single antenna 
or two antennas each to give ur = 3 or ; = 6 respectively. 
In the shaded coverage area of this edge-excited cell, we drop 
three or four users uniformly in space giving N = 3 or 
N = 4. For each user, lognormal shadow fading and path loss 
is considered, where the standard deviation of the shadowing 
is 8dB and the path loss exponent is 7 = 3.5. The transmit 
power of the sources is scaled so that the best signal received 
at the three BS locations is greater than 3dB at least 95% of 
the time. Even though the analysis in this paper is valid for 
any set of channel powers, the above methodology allows us 
to investigate the accuracy of the performance matrices for 
realistic sets of channel powers. 

In Figs. [3] H] and [5] the case of three single antenna users 



and three distributed BSs with a single receiver antenna is 
considered. Here, we investigate both the approximate SINR 
distributions and the approximate SER results for an MMSE 
receiver. In Fig. [3] the approximate cdfs of the output SINR 
are plotted alongside the simulated cdfs. Results are shown for 
four random drops and, the results are for a particular user (the 
first of the three). The agreement between the cdfs is shown 
to be excellent. Note that this agreement is good across all 
drops, from Dl which has a very poor SINR performance 
to D4 with a much higher SINR performance. The use of 
physically motivated drops rather than ad-hoc scenarios is 
useful as it assesses the accuracy of the analysis in plausible 
channel conditions. 

In Fig. |U the approximate SER curve is plotted alongside the 
simulated values. Results are shown for three drops and QPSK 
modulation. The agreement between the SER results is shown 
to be excellent across all three drops at SERs below 10~ 2 . 
Again, this agreement is observed over a wide range with Dl 
having much higher SERs than D3. In Fig.|4]and also in Figs. 
[6]|7] the SER is plotted against the transmit SNR, 7. This is 
chosen instead of the receive SNR to separate the curves so 
that the drops are visible and are not all superimposed, which 
tends to happen when SER is plotted against receive SNR. 
In Fig. [5] the approximate cdfs of the SNR are plotted 
alongside the simulated cdfs for a ZF receiver. Results are 
shown for four random drops. This is the companion plot to 
Fig. [3] with the same system but a ZF receiver rather than an 
MMSE receiver. The accuracy of the results in Fig. [3] and Fig. 
[5] is interesting, especially when you observe that the Fig. [3] 
analysis uses d7Tl i. a simple mixture of 3 exponentials, and 
Fig. uses d43l which is a single exponential in this case. 
In Fig. [6] and [7] the case of four single antenna users and 
six distributed receive antennas (two at each BS location) is 
considered. High SNR SER curves are plotted alongside the 
simulated values. Results are shown for both MMSE (Fig. 
[6]l and ZF (Fig. IT) with QPSK modulation. The agreement 
between the simulated SER and the high SNR approximation 
is shown to be less accurate than in Fig. |4] with very close 
agreement requiring low error rates around 10 ~ 4 . This is un- 
surprising, as the greater number of system dimensions gives 
greater freedom for the channel powers to vary substantially 
over the links. 



VIII. Conclusion 

The performance of MMSE and ZF receivers in a macrodi- 
versity layout is a long-standing, unsolved research problem. 
In this paper, we make the first major progress towards solving 
this problem for the general case of an arbitrary number 
of transmit and receive antennas. The analysis is based on 
a derivation which targets the characteristic function of the 
output SINR. This leads to an expected value which is highly 
complex in its exact form, but can be simplified by the use of 
an extended Laplace type approximation. This methodology 
is able to produce approximate results for both the SINR 
distribution and the SER which are simple enough to be 
employed as the basis of a scheduling algorithm. The SINR 
distribution is shown to have a remarkably simple form as 
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Fig. 2. Network MIMO/edge-excited cell scenario where three base stations 
serve users in a three-sector cluster. To reduce clutter, only two users are 
shown. 




SINR [dB] 



Fig. 3. Approximate and simulated SINR cdf results for the TV = 3, tir = 3 
scenario. Results are shown for the first of three users for four arbitrary drops 
and a MMSE receiver. 

a generalized mixture of exponentials. Also, the asymptotic 
SER results produce a remarkably compact metric which 
captures a large part of the functional relationship between 
the macrodiversity power profile and SER. In addition, since 
both SINR and SER results are available, scheduling can be 
performed on the basis of either metric. 

Appendix A 
Calculation of \D\ Perm (p~ 1 Q 2 s J 

The permanent of the denominator in (130b can be expanded 

as 

Perm (d^Q^j = ^perm (jo^Q^ ^ ^ , (92) 

where (JN-i,n a is an ordered subset of {71r} = 
{1,2,..., nji} of length N — 1 and the sum is over all (>£f -J 
such subsets. Noting the fact that perm (EX) = |E perm (X), 




5 10 15 20 25 30 
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Fig. 4. Approximate and simulated SER results for the TV = 3, tir = 3 
scenario with QPSK modulation. Results are shown for the first of the three 
users for three arbitrary drops and a MMSE receiver. 




SNR [dB] 



Fig. 5. Approximate and simulated SNR cdf results for the TV = 3, ur = 3 
scenario. Results are shown for the first of three users for four arbitrary drops 
and a ZF receiver. 



for a square diagonal matrix E and ( TTTI i. d92l can be further 
simplified to give 

perinffQ,)'"- 11 ) 

U <XN-l, nft 

Using d93l l, the denominator in ( f30b becomes 
|D| Perm (ir 1 ^) =Y\ D ^ B 

xperm((Q 2 )^: i 1 >J, (94) 

where o"l.ur is the ordered subset of length L of 
{1,2,..., ur] which does not belong to <TN-i,n R an d L = 
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Calculation of E ■ 



oH + H^H 2 



-* — Simulation 
» High SNR approx 



Fig. 6. Approximate and simulated SER results for N = 4, tir = 6, i.e., 
two receive antennas at each BS with QPSK modulation. Results are shown 
for the first of four users for three arbitrary drops and a MMSE receiver. 




-• — Simulation 

■* High SNR approx. 



Fig. 7. Approximate and simulated SER results for the N = 4, ur = 6, 
i.e., two receive antenna at each BS scenario with QPSK modulation. Results 
are shown for the first of four users for four arbitrary drops and a ZF receiver. 



Ur — N + 1. Expanding \D Sl 

L '~3t 



gives 



Substituting d95t in d94b gives the desired result 



Pi 



where 

Pi = 



\D\Vecm(p- 1 Q 9 ) = ^2(-jty 

t=0 

E T ^((^) ff£ ,„Jpenn((Q 2 )^J(a 2 



(95) 



(96) 



(97) 



Similar expectation results for random determinants can 
also be found in lfl4l . However, for completeness we present 
the particular result needed for the MMSE analysis here. Let 
Ai, A2, . . • , Ajv-i be the ordered eigenvalues of H 2 H 2 . Since 
n>R > N, all eigenvalues are non zero. Then 



E- 



a l I + H" H 2 



}=s|n( ( x 2 + A 4 )} 



(98) 



'jV-1 



2 \JV-t-l| 



i=0 



(99) 



where (|99l is from (fT~3T > and Lemma [2] Therefore, the building 
block of this expectation is E I Tr^ [H 2 H% ) \ ■ From Lemma 

m 



h 2 h 2 



Tr 4 (iff H 2 ) = 

a 

Therefore, from Lemma \T\ 

Efa (fTfiTa)} =^Perm((g 2 r ~- 1 ), 

(7 

where the ur X (N — 1) matrix, Q 2 , is given by 

E{H 2 oH 2 }=Q 2 . 
Then, the final expression becomes 

E { \a 2 I + H §H 2 J } =J2Y, Perm ((Q2) Cri ' N - 1 ) 

i=0 <r 



(100) 



(101) 



(102) 



N-i-l 



(103) 



Appendix C 
Calculation of \D\ E j cr 2 / + H"D 1 H 2 1 1 

As a simple extension of expectation in the numerator of 
J, the expectation in the denominator can also be calculated 
to get 

N-l 

E^I + H^D-'H^^^i-jt) (a 2 )™, 

k=0 

(104) 

where 



Vfc(-it) = X) Perm (( £r1 ^ 



and 



^0 (-J<) = 1. 



(105) 



(106) 



The term in d 1 05t can be simplified using (jT7j and Corollary 
|2] to obtain 



4>k (-jt) = ■ 



Perm 



(«? 2 C;; } 



D 



(107) 



Cfc.n 
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Then, 



N-l 



\D\ E{\a*I + H^D- l H 2 \)=Y, & (-#) 



fc=0 



(108) 



where ^ (— it) = |JJ| ^fc (~jt). From (11071 ). we can get 



& (-it) = £ |fl VM , Perm ((Q 2 )^ 1} 



(109) 



where ay, 



is a length n# — k subset of {1, . . . ,n,f>} 



which does not belong to <7k,n R - Therefore, it is apparent that 
£fc (~jt) is a polynomial of degree ur — k. It is clear from 
(ITOBl that, when cr 2 = 0, (TT081 collapses to ([94}. Clearly, 
\D\ E | cr 2 J + H^D ~ l H% I is a polynomial of degree ur, 
as £0 ( — it) = l-D I is me highest degree polynomial term in 



(U081 Then, 
Hence 



riR — k 



= £(^)'M™**->J- (no) 



rin — k 



(7 i=Q 



x^nn(«3 2 )<»-'>), 



so that (—it) becomes 



riR-k 



& (-it) = £ ( "4" I ^ 



i=0 

E 

i=0 



(HI) 

(112) 
(113) 



where 





* = E Tr * (( p ik R -*, J Perm > (114) 



and simplifies to give 

<Pi0 = Tr 4 (Pi) . 
Equation ( II 13b follows from the fact that 

T^((Pi) ffnfl _ feR ) =0 for i>n R -A;. 



(115) 



(116) 



Therefore, (11081 l can be written as 



JV-l rij 



} = E E^ - - 7 '*)*^ 



fc=0 i=0 



/ o\N—i—k—l 

x (a 2 ) 



(117) 



which is in turn can be given as 



\d\e{ 



a z I + H"D 



} = E(-i*)Vi> (ii 8 ) 



where 



iV-l 



= & k i a ' 



.N-i-k-l 



(119) 



fc=0 



i=0 
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